Polarity-induced oxygen vacancies at LaA103|SrTi03 interfaces 
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Using first-principles density functional theory calculations, we find a strong position and thickness 
dependence of the formation energy of oxygen vacancies in LaA103|SrTi03 (LAO|STO) multilayers 
and interpret this with an analytical capacitor model. Oxygen vacancies are preferentially formed at 
p-type SrO|AlC>2 rather than at n-type LaO|Ti02 interfaces; the excess electrons introduced by the 
oxygen vacancies reduce their energy by moving to the n-type interface. This asymmetric behavior 
makes an important contribution to the conducting (insulating) nature of n-type (p-type) interfaces 
while providing a natural explanation for the failure to detect evidence for the polar catastrophe in 
the form of core level shifts. 



PACS numbers: 68.35.-p, 68.35. Ct, 73.20.-r 

Introduction: Extremely high carrier mobilities have 
recently been observed when interfaces consisting of LaO 
and TiC>2 layers are formed between insulating LaA103 
and SrTiOa perovskites [1] . Even though the physical ori- 
gin of this metallic behavior is still under debate [2-7], 
most experimental [8-19] and theoretical [20-28] studies 
have reached a consensus that the so-called polarity dis- 
continuity between these materials plays a crucial role; in 
the absence of any relaxation mechanism, alternate stack- 
ing of positively (LaO + ) and negatively (A10^~) charged 
layers on the non-polar STO substrate would give rise to 
a huge effective internal electric field, leading to a diver- 
gence of the electrostatic potential with increasing thick- 
ness of LAO. Three mechanisms have been suggested to 
avoid this instability: charge transfer [2], atomic relax- 
ation [12, 19-23], or the creation of oxygen vacancies [4, 5] 
and other defects [7, 13]. 

The first mechanism refers [2] to the transfer of elec- 
trons from a surface AIO2 layer to the interface TiC>2 
layer by the internal electric field. The excess charge at 
the interface balances the polar discontinuity and leads 
to conducting behavior of the interface. This mechanism 
is strongly supported by the observation of an insulator- 
metal transition induced by either an external electric 
field or by increasing the thickness of the LaAlC>3 layer 
[9]. However, direct experimental evidence of charge 
transfer, in the form of core level shifts, has not yet been 
found [6] . The insulating behavior of the p-type interface 
[1] is also not readily accommodated in this picture. 

The second mechanism, atomic relaxation in the pres- 
ence of the internal electric field, that is analogous to the 
buckling of Ti-O-Ti chains in an external field in SrTiC>3, 
has been discussed by a number of authors [12, 19-23]. 
It can eliminate the diverging potential by introducing 
a compensating electric field. A third way to resolve 
the polar instability is to introduce defects at interfaces. 
Oxygen vacancies (and other defects) created during the 
growth of LAO on STO arc invoked by Herranz et al. [4] 
and Kalabukhov et al. [5] to understand the high mo- 



bility carriers. The long-relaxation time of the electric- 
field-induced insulator-metal transition [15] suggests the 
possibility of interface defect diffusion. 

While the first two mechanisms have received much 
theoretical attention [20-27], the relationship between 
the creation of defects and polarity has not been clarified 
theoretically. To demonstrate the coupling of polarity 
and oxygen vacancy formation and throw some light on 
the interplay with atomic relaxation and charge transfer, 
we calculate from first-principles the formation energy of 
oxygen vacancies in LAO|STO multilayers as a function 
of their location in the multilayer. 

Method: We focus on (m,m) LAO | STO multilayers 
containing m layers each of LAO and STO with alternat- 
ing p- and n-type interfaces. Because samples are grown 
on STO substrates, we fix the in-plane lattice constant 
at the calculated equilibrium value of STO and calculate 
the out-of-plane lattice constant by minimizing the total 
energy of strained bulk LAO. Oxygen vacancies are mod- 
elled in a 2 x 2 lateral supercell and for each vacancy posi- 
tion all atoms are allowed to relax with the volume of the 
structure fixed. Most of the results reported below were 
obtained with the 159 atom supercell depicted in Fig. 1 
containing a p and n interface. The periodically repeated 
single oxygen vacancy in a layer consisting of 2 x 2 unit 
cells should be compared to the ~25% oxygen vacancy 
concentration in a layer suggested by experiment [2] . The 
local density approximation (LDA) calculations were car- 
ried out with the projector augmented wave method [29] 
as implemented in the Vienna Ab-initio Simulation Pack- 
age (VASP) [30]. A kinetic energy cutoff of 500 eV was 
used and the Brillouin zone of the 159 atom supercell was 
sampled with an 8 x 8 x 2 k-point grid in combination 
with the tetrahedron method. 

The concentration of vacancies in thermodynamic 
equilibrium is determined by the vacancy formation en- 
ergy, Eq^. Since we want to focus on how it depends on 
the position (z) of the defect for a charge neutral system, 
we express Eq^ c for an oxygen vacancy in an LAO | STO 
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FIG. 2: Position dependence of the formation energy of an 
oxygen vacancy in a 2 x 2 x (4 + 4) LAO|STO multilayer cal- 
culated from first principles without relaxation (symbols) and 
using an analytical capacitor model (solid line), n and p in- 
terfaces are indicated by vertical red and blue dashed lines. A 
schematic diagram of the capacitor model is shown in the up- 
per panel. The electrostatic potential profile for the vacancy- 
free structure is shown as a dotted line. The vertical black 
line at a distance d from the n interface represents the oxy- 
gen vacancy layer. Two excess electrons are transferred to the 
TiC>2 layer at the n interface (shaded grey line). 



FIG. 1: Left panel: the unit cell of a 2 x 2 x (4 + 4) LAO|STO 
multilayer with an oxygen vacancy at the p-type interface. 
Blue spheres represent oxygen atoms and the oxygen vacancy 
is marked by a white sphere. Charge density isosurfaces cor- 
responding to a value of 0.015 e/A 3 for occupied states in 
the conduction bands are coloured red. Right panel: plane- 
averaged charge density as a function of z for oxygen vacancies 
at p- (red) or n- (blue) type interfaces. 

multilayer as the difference in total energies of super- 
cells with and without a vacancy, Eq^ c = EYq c — -Esc 
[31]. Because it is useful to relate this to the forma- 
tion energy of an oxygen vacancy in bulk STO, E^q, 
all of the vacancy formation energies discussed below are 
AE%* C = E^ c - Ej%° . Though E^ C Q converges only 
slowly with the size of supercell [32], this does not al- 
ter the conclusions we will draw; E^q c spans a range of 
±1.0 eV so the effect of polarity is much larger than the 
uncertainty in E^q which is converged to ±0.09 eV. 

Unrelaxed results: We begin by calculating AEq^(z) 
without atomic relaxation. The most striking feature of 
the results shown in the lower panel of Fig. 2 is the asym- 
metry for forming a vacancy at the n and p interfaces. 
A£qJ c is lowest when the oxygen vacancy is at the p in- 
terface, highest close to the n interface and is nonlinear 
in z. It spans a range of about 3 eV between the two in- 
terfaces and differs in the LAO and STO layers. We can 
capture the essential behaviour of AEq^(z) in terms of 
a modified parallel-plate capacitor model. 



Model: The average electrostatic potential of defect- 
free LAO | STO multilayers, as probed by the energy lev- 
els of core states in LDA calculations, exhibits a simple 
symmetric triangular form as if the n and p interfaces 
were positively and negatively charged with charge den- 
sity ±<7 = Te/a 2 where a is the lattice constant of bulk 
STO. As sketched in the top panel of Fig. 2, the plates 
of the capacitor are separated by a thickness d\ (d 2 ) of 
insulating LAO (STO) with dielectric constants £\ (£2) 
determined by the electronic polarization only in the ab- 
sence of ionic relaxation. Such a model was recently used 
to describe the evolution of the dielectric properties of 
LAO I STO multilayers with increasing layer thickness re- 
sulting in an insulator-metal transition [26]. Based on 
this simple capacitor model, the estimated internal elec- 
tric fields are huge, |l.2 x l0 11 V/m or ~ 0.9 V/unit cell 
in LAO, and the electrostatic potential (dotted line in 
Fig. 2) diverges with increasing thickness of LAO. 

We extend this model to encompass the layer of va- 
cancies constructed in our supercell approach. Because 
oxygen is divalent, removal of a neutral oxygen atom in a 
bulk insulating material such as LAO or STO leaves two 
excess electrons in the conduction band weakly bound 
to the oxygen vacancy. In an LAO | STO multilayer, the 
potential energy of the electrons in the internal electric 
field far exceeds this binding energy and the total en- 
ergy can be reduced by moving the two electrons to the 
conduction band minimum at the n interface leaving a 
sheet of positive charge at the oxygen vacancy plane a 
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FIG. 3: Position dependence of the formation energy of an 
oxygen vacancy in a 2 x 2 x (4 + 4) LAO|STO multilayer 
relative to that of a vacancy in bulk STO (horizontal dotted 
line) calculated from first principles with relaxation. Black 
and grey symbols are for vacancies in AO and BO2 layers 
respectively. The formation energy of an oxygen vacancy in 
bulk LAO is shown as a dashed horizontal line. 



distance d away. Assuming that the excess electrons are 
on Ti ions at the Ti02|LaO interface, independent of 
where the oxygen vacancies were formed, we can calcu- 
late the d dependent change in the electrostatic energy to 

be d 1 J 2 +d° El {— ^1 + fj(^2 — d)}d. It comprises two parts: 
the energy to insert a sheet of positive charge density a/2 

2a 2 le 

in the LAO|STO capacitor background, — die9+rf ° ei did, 
and the potential energy of the positively and negatively 

charged sheets d °+d e {d\ + §J(da ~ d)}d. The calcu- 
lated LDA core level shifts of ~ 0.9V/unit cell can be 
used to estimate e\ + £2 ~ 52 leaving one free parameter 
in the model, the ratio £2/^1 ■ Taking this to be 1.5 results 
in the solid curve in Fig. 2. The good fit of this simple 
model makes it clear that the internal fields induced by 
the polar layered structure can lower the formation en- 
ergy of oxygen vacancies at the p interface very substan- 
tially and that the origin of the asymmetry in formation 
energies in LAO and STO is the difference in their di- 
electric constants. The residual interaction between the 
field-ionized oxygen vacancies and electrons accounts for 
the nonlinear behavior of the formation energy. 

Atomic relaxation can be expected to strongly sup- 
press the polarity. Nevertheless, when our structures are 
fully relaxed, some essential features of AEq^(z) are un- 
changed, see Fig. 3. In particular, the formation energy 
has a minimum at the p interface and a maximum close 
to the n interface while the minimum formation energy 
is more than 1 eV lower than in bulk STO. Including 
atomic relaxation differentiates between vacancy forma- 
tion in AO and BO2 layers; the latter are energetically 
more favourable though the behaviour as a function of 
z is essentially the same. The lower formation energy 
in BO2 layers can be understood in terms of the types 
of relaxations possible within the constraints imposed by 
stacking the different layers in a multilayer. To simplify 
the discussion, we focus on the less favourable case of 
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FIG. 4: Upper panel: projection of the Ti-O-Ti separation 
along the z direction due to buckling for clean interfaces (A) 
and when vacancies are formed at the p (□) and at the n 
(O) interface. Lower panel: formation energy of an oxygen 
vacancy as a function of the multilayer thickness for p and 
n interfaces in a 2 x 2 x (m + m) LAO | STO multilayer with 
relaxed structure. 



oxygen vacancies in AO layers so that our conclusions 
will also be applicable for B02-layer vacancies. 

Critical thickness: For a capacitor with fixed charge 
density ex, the electric field is constant and the electro- 
static potential increases as the plate separation (LAO 
thickness) is increased, a feature that is supported by 
both experimental [9, 15] and theoretical [20-23, 26] stud- 
ies. Because of its dependence on the electrostatic poten- 
tial, we expect the formation energy of oxygen vacancies 
to depend on the multilayer thickness. Since the mini- 
mum and maximum formation energies occur at or close 
to the p and n interfaces, respectively, we focus on these 
formation energies. We further assume equal thicknesses 
m of the STO and LAO layers and plot AEq 1 ^ as a func- 
tion of m in Fig. 4. At the n interface it is almost constant 
in value and ~ 0.6 eV higher than for bulk STO. At the 
p interface however, it decreases with increasing m and 
becomes negative for a critical value of m between 3 and 
4. Thus, for m > 4 oxygen vacancies are preferentially 
formed at p interfaces rather than in the bulk of the ma- 
terials or at the n interface. To better understand this 
behaviour, it is useful to analyze how the atomic relax- 
ation depends on the polarity of the system. 

In the absence of atomic relaxation, uncompensated in- 
ternal electric fields in an LAO | STO multilayer (as mea- 
sured for example by core level shifts) lead to the occu- 
pied LAO oxygen p band rising in energy from a mini- 
mum at the n interface to cross the unoccupied Ti d con- 
duction band minimum so the interface becomes conduct- 
ing. When atomic relaxation is allowed, so-called zigzag 
buckling of the Ti-O-Ti and Al-O-Al chains in the Ti0 2 
and AIO2 planes occurs [12, 20-23, 26] that compensates 
the electric fields and makes the whole system insulat- 
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ing again. However this buckling itself costs energy and 
as m increases, oxygen vacancy formation may become 
a more favourable way of compensating the polarity and 
releasing the strain. For p-n multilayers, vacancies at the 
p interface will suppress the internal electric field much 
more strongly than at the n interface, involve much less 
Ti-O-Ti buckling and be energetically favourable for m 
larger than some threshold value. This argument is sup- 
ported by a plot of the average relative displacements 
(excluding interfacial layers) of the Ti ions and O ions in 
the z direction shown in Fig. 4. The buckling is clearly 
strongly suppressed when oxygen vacancies are formed at 
the p interface and reaches a saturation value when the 
multilayer thickness is increased. However, compared to 
the clean interface case, it is almost unaffected by in- 
troducing oxygen vacancies at the n interface. The core 
level shifts are reduced from 0.9 to 0.15 V/unit cell. 

Electronic structure: A single oxygen vacancy donates 
two electrons to the system and the distribution and char- 
acter of these excess electrons will dominate its transport 
properties. In a clean LAO|STO multilayer, the electro- 
static potential at the p interface is lower than at the 
n interface so even if oxygen vacancies are generated at 
the p interface, the excess electrons will be driven by the 
electrostatic potential to the n interface. The charge dis- 
tribution of the occupied conduction band states is plot- 
ted in Fig. 1. Consistent with previous studies [17, 24], 
these have Ti d^-orbital character if the Ti ions are close 
to the interface, otherwise they consist of a mixture of 
Ti d xz - and dj, z -orbital character. Though the introduc- 
tion of oxygen vacancies represents a major change to the 
atomic structure, the change to the electronic structure 
is minor compared to a clean LAO|STO multilayer with 
only ro-type interfaces. 

Even though it is energetically favourable to form oxy- 
gen vacancies at a p interface, the excess electrons are 
transferred to the n interface which will tend to be con- 
ducting while the p interfaces will be insulating. This 
asymmetry implies a spatial separation between impurity 
scattering and transport regions which, by analogy with 
proximity doping in semiconductor heterostructures, can 
give rise to a high mobility at LAO|STO interfaces. This 
result agrees with recent experimental data [2] which 
shows evidence for the presence of oxygen vacancies at 
the p interface while the mobile carriers are at the n in- 
terface. Though most experimental studies have been 
made [1, 2, 9] on samples consisting of STO substrates 
covered with several layers of LAO which only contain a 
single n interface, the conclusions of this paper are quali- 
tatively valid for the single interface case if we regard the 
surface of LAO as a pseudo-p type interface. 

Conclusion: Using first-principles calculations, we 
show how oxygen vacancy formation, charge transfer and 
atomic relaxation in response to polar discontinuity at 
LAO | STO interfaces are strongly coupled. Oxygen va- 
cancies are preferentially formed at p rather than n inter- 



faces and the thickness-dependent formation energy pro- 
vides an alternative explanation for the critical thickness 
observed in experiments while simultaneously explaining 
the failure to observe core level shifts. The conduction 
electrons produced when an oxygen vacancy is formed 
at a p interface move to the n interface where their in- 
teraction with the vacancies is minimal explaining the 
observed high mobilities. 
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